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I.  INTRODUCTION 


Contract  N00014-82-K-0805  between  Applied  Research  Laboratories, 
The  University  of  Texas  at  Austin  (ARL:UT),  and  the  Office  of  Naval 
Research  (ONR)  provided  for  research  on  nonlinear  effects  in  long  range 
underwater  propagation.  The  period  of  the  contract  was  1  September 

1982  -  30  November  1984. 

The  general  purpose  of  the  research  was  to  determine  the  extent  to 
which  nonlinear  effects  are  important  in  long  range  underwater 
propagation.  Pressure  pulses  of  the  sort  produced  by  explosives  were 

considered.  In  order  to  work  on  several  aspects  of  the  problem  at  the 

same  time,  we  divided  the  research  into  three  tasks. 

Task  I.  Shock  pulse  propagation  in  a  homogenous  ocean.  In  this 
task  the  effect  of  inhomogeneity  of  the  ocean  is  suppressed,  and 
the  pulse  is  treated  simply  as  a  spherically  spreading  wave. 
Ordinary  nonlinear  effects--the  decay  of  the  shock  at  the  head  of 
the  pulse  and  the  waveform  distortion  of  the  tail  behind--are 
included  as  are  real-ocean  attenuation  and  dispersion. 

Task  II.  Nonlinear  propagation  in  a  depth  dependent  ocean.  Here 
the  inhomogeneity  of  the  ocean,  in  particular,  stratification,  is 

taken  into  account.  The  primary  concern  is  the  effect  of 
inhomogeneity  on  finite-amplitude  distortion.  Shock  effects  per  se 
and  ordinary  attenuation  and  dispersion  are  ignored. 

Task  III.  Nonlinear  propagation  in  a  caustic  region.  In  this  task 
the  special  problem  of  the  passage  of  a  finite-amplitude  pulse 
through  a  caustic  region  is  taken  up. 

In  the  end,  of  course,  it  would  be  desirable  to  combine  the  results  of 
the  three  tasks.  We  wished  to  be  able  to  calculate  the  propagation  of  a 
finite-amplitude  pulse  through  a  real,  inhomogeneous  ocean,  even  when 
the  propagation  path  includes  caustics. 


For  administrative  purposes,  a  small  additional  project  on 
nonlinear  mixing  of  two  beams  of  sound  was  later  added  to  the  contract. 

Although  D.  T.  Blackstock  was  the  principal  investigator  for  the 
contract,  the  primary  senior  person  who  planned  and  carried  out  the 

research  was  C.  L.  Morfey.  A  faculty  member  at  Institute  of  Sound  and 
Vibration  Research  (ISVR),  University  of  Southampton,  Morfey  did  most  of 
his  work  while  on  leave  at  ARL:UT  during  the  first  half  of  1983.  At 
other  times  during  the  life  of  the  contract,  he  pursued  the  research  as 
a  consultant.  R.  Buckley,  also  a  consultant  from  ISVR,  performed  the 
work  on  the  problem  of  caustics  (Task  III).  ARl :UT  staff  members  who 
assisted  Morfey  were  J.  M.  Estes  (development  of  numerical  propagation 
programs  for  Task  I),  T.  L.  Foreman  (advice  on  the  ray  tracing  program 
for  Task  II),  and  J.  Hardin  and  D.  Gring  (computation  assistance).  The 
separate  project  on  nonlinear  mixing  of  sound  beams  was  the 

responsibil ity  of  consultant  P.  J.  Westervelt  (Brown  University). 

Scientific  Officers  for  ONR,  Code  425  Underwater  Acoustics,  were, 
successively,  P.  H.  Rogers,  W.  I.  Roderick,  and  R.  M.  Fitzgerald. 


II.  TECHNICAL  DISCUSSION 


An  enormous  amount  of  literature  exists  on  long  range  propagation 
of  sound  in  the  ocean.  Much  broader  than  it  is  deep,  the  ocean  is  a 
natural  waveguide,  albeit  a  very  complicated  one.  The  medium  is 
inhomogeneous  and  the  boundaries  are,  to  say  the  least,  nonideal.  The 
three  main  techniques  that  have  been  used  in  analytical  studies  of  long 
range  propagation  are  ray  theory,  the  parabolic  approximation,  and 
normal  mode  theory.  All  are  based  on  the  linear  wave  equation. 
Although  nonlinear  effects  are  rarely  considered,  it  is  interesting  to 
note  that  two  of  the  techniques,  ray  theory  and  the  parabolic 
approximation,  may  be,  and  in  fact  have  been,  generalized  to  include 
nonlinear  propagation  effects.  It  would  be  exceedingly  difficult  in  any 
general  way  to  modify  normal  mode  theory  to  Include  nonlinearity. 


Except  for  studies  of  finite-amplitude  distortion  of  the  shock 

pulse  produced  by  an  underwater  explosion, 1-4  little  interest  has  been 

shown  in  possible  effects  of  nonlinearity  on  long  range  propagation. 

The  prevailing  opinion  seems  to  be  that,  while  perhaps  important  near 

the  source,  where  the  amplitude  is  high,  nonlinear  effects  surely  become 

insignificant  for  ranges  beyond  a  few  hundred  meters  or  so.  This  belief 

3 

persists  despite  the  fact  that,  as  Marsh,  Mellen,  and  Konrad  pointed 

2 

out  in  1965,  Arons'  amplitude  scaling  law  is  known  to  hold  over  five  or 
six  decades  in  range.  This  law  is  that  the  peak  amplitude  ps  of  a  wave 
from  an  explosion  decays  with  range  R  as 

p$  ®  R'1*13  (1 

(not,  as  given  by  linear  theory,  as  ps«R-*).  A  second  and  more  telling 
argument  against  considering  nonlinear  effects  has  to  do  with  frequency. 
Because  attenuation  increases  rapidly  with  frequency,  many  long  range 
measurements  are  limited  to  low  frequencies,  and  low  frequency  signals 
are  known  to  be  highly  resistant  to  finite-amplitude  distortion.  Two 
counter  arguments  are  as  follows:  (1)  nonlinear  distortion  is 

cumulative  and,  though  small  over  short  distances,  the  distortion 
effects  may  add  up  to  produce  appreciable  effects  over  long  ranges; 


V"'  »'»«■*  k  ^  . 


(2)  studies  of  both  finite-amplitude  tones  and  finite-amplitude  noise 
have  shown  that  for  the  high  frequency  end  of  the  spectrum  no  range 
exists  at  which  linear  theory  is  adequate.5  As  an  example  of  (2), 
consider  a  spherically  spreading  sound  field  generated  by  a  pure  tone 
source.  It  has  been  found  that  the  amplitude  of  the  nth  harmonic 


distortion  component  pn  varies  as 

_  -na.R 
p„  -  R-"  e  1 

not,  as  expected  from  linear  theory, 

-i  -ct_  R 
Pn  «  R  e  n 


(2) 


(3) 


Here  and  otn  are  the  attenuation  coefficients  at  the  frequencies  of 
the  fundamental  and  the  nth  harmonic,  respectively. 


About  a  year  before  our  work  on  the  problem  of  finite-amplitude 
propagation  in  the  ocean  started,  B.  E.  McDonald  and  W.  A.  Kuperman  at 
Naval  Ocean  Research  and  Development  Activity  (NORDA)  began  to  attack 
the  problem  by  using  the  (time-domain)  nonlinear  parabolic  equation.5 
This  equation  is  well  suited  to  numerical  solution,  and  their  work  may 
be  characterized  as  computationally  intensive.7"*0  Taken  into  account 
are  nonlinear  propagation  distortion,  refraction  (including  caustics), 
and  diffraction  (they  considered  a  beam  of  sound,  as  is  appropriate  when 
the  parabolic  equation  is  used).  Not  taken  into  account  is  ordinary 
absorption  by  the  medium;  however,  shock  losses  are  included.  Among 
their  reported  results  are  the  evolution  of  waveforms  and  the  behavior 
at  caustics,  in  particular  the  effect  of  nonlinearity  on  phenomena  at 
caustics . 


We  decided  to  use  ray  theory  for  our  investigation.  Ray  theory  is 
simple  and  versatile  and  is  appealing  for  physical  interpretation.  The 
extension  to  nonlinear  acoustics  is  relatively  well  known.  Finally,  the 
analytical  work  can  be  carried  a  great  deal  further  with  ray  theory  than 
with  the  parabolic  equation.  Although  numerical  computation  must  in  the 
end  be  used  to  calculate  waveforms,  the  computation  is  in  principle  very 
simple.  Moreover,  the  use  of  ray  theory  allows  one  to  easily 
investigate  the  individual  effects  of  the  various  physical  mechanisms. 


such  as  nonlinearity,  absorption,  and  refraction,  on  the  propagation. 
The  only  significant  disadvantage  associated  with  ray  theory  is  that 
caustics  pose  a  difficult  problem.  As  already  noted,  we  planned  to  deal 
with  caustics  as  a  separate  task. 

The  two  approaches,  ours  and  that  of  the  NORDA  group,  were 
compl ementary  and  provided  complementary  results. 


A.  Task  I.  Shock  Pulse  Propagation  in  a  Homogenous  Ocean 

In  this  task  nonlinear  propagation  distortion  (described  by  weak 
shock  theory)  and  real-ocean  attenuation  and  dispersion  are  considered 
(inhomogeneity  of  the  medium  is  considered  in  Task  II,  not  here).  Shock 
pulse  parameters  of  particular  interest  are  the  peak  pressure  ps  and  the 
decay  time  behind  the  shock  fls. 

The  objective  of  Task  I  was  to  construct  a  numerical  model  of  long 
range  shock-pulse  propagation  which  would  allow  p$  and  6S  to  be 
calculated  as  a  function  of  range  R.  The  calculation  would  start  with 
an  assumed  initial  waveform.  Our  intention  was  to  compare  the  numerical 
model,  in  which  relaxation  attenuation  and  dispersion  are  included  along 
with  nonlinear  waveform  steepening,  with 

(1)  experimental  data  as  summarized  in  the  Arons  scaling  laws,  and 

(2)  weak-shock  theory,  which  does  not  include  relaxation  effects. 

Particular  questions  of  interest  were  the  following. 

•  Why  do  the  empirical  scaling  laws  for  ps  and  G$  fit  the 
available  data  better  than  an  exact  solution  based  on 
weak-shock  theory? 

•  Out  to  what  distance  may  weak-shock  theory  be  used  to  describe 
the  shock  pulse? 

•  What  influence  does  the  first  bubble  pulse  have  on  the  shock 
pulse  parameters  at  long  range?  (According  to  the  weak-shock 
model,  there  is  no  interaction.) 
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•  Beyond  what  distance  may  nonlinear  effects  be  neglected,  for 
purposes  of  describing  the  evolution  with  range  of  certain 
waveform  parameters? 

The  propagation  algorithm  that  was  developed  is  closely  related  to 
one  devised  by  Pestorius^  for  airborne  propagation  of  finite-amplitude 
waves  through  a  tube.  Modifications  were  made  to  allow  for  spherical 
spreading  and  a  switch  of  the  acoustical  medium  to  sea  water.  One 
starts  with  the  time  waveform  of  the  wave  at  the  source  and  computes  the 
change  in  waveform  as  incremental  steps  in  range  are  taken  (when  the 
wave  spreads  spherically  the  steps  are  equal  increments  in  In  R).  For 
each  step  the  wave  is  first  subjected  to  ordinary  attenuation  and 
dispersion  (appropriate  for  ocean  conditions)  over  that  step.  This 
operation  is  carried  out  in  the  frequency  domain.  After  a 
transformation  back  to  the  time  domain  is  made,  the  waveform  is 
corrected  for  nonlinear  propagation  distortion  by  using  the  relations  of 
weak  shock  theory.  The  sequence  of  absorption  and  dispersion  followed 
by  nonlinear  distortion  is  repeated  for  each  step  until  the  desired 
propagation  distance  is  reached. 

The  algorithm  was  used  to  perform  numerical  propagation 
"experiments"  for  equivalent  TNT  charge  weights  in  the  range  0.6  g  to 
20  kg.  The  following  initial  waveforms  were  used. 

(1)  P  =  e~T  (simple  exponential  tail),  (4) 

(2)  P  =  [l-(l-k)T]e'kT  (Brode12),  (5) 

1 3 

(3)  A  modified  version  of  Wakeley's  waveform  with  the  first 
bubble  pulse  included. 

Here  P=Rp/RQpm,  T=(t-R/c)/0m,  p  is  acoustic  pressure,  R0>  p  ,  and  9m  are 
the  starting  values  of  R,  ps,  and  0  ,  respectively,  t  is  time,  and  c  is 
sound  speed;  Brode's  constant  k  is  3/3. 


Rogers  used  weak-shock  theory  to  obtain  an  analytical  solution  for 
waveform  (1).  Ordinary  absorption  and  dispersion  were  not  included.  In 
a  sense  our  effort  under  Task  I  is  an  extension  of  Rogers’s  work,  since 
we  considered  a  more  realistic  ocean  and  more  realistic  waveforms. 
Indeed  we  were  able  to  explain  some,  but  not  all,  of  the  differences 
Rogers  noticed  when  he  compared  his  analytical  solution,  Arons's  scaling 
law,  and  measured  data.  For  details,  see  References  84-2  and  85-1.* 

The  abstract  from  85-1  summarizes  the  results  under  Task  I. 

"The  propagation  of  idealized  model  waveforms,  representative  of 
underwater  explosions  measured  at  moderate  distances  (particle  Mach 
number  =  0.06),  has  been  simulated  numerically  out  to  large 
distances  in  an  idealized  uniform  ocean.  Besides  chemical 
relaxation,  the  calculation  procedure  includes  nonlinear  waveform 
distortion  and  thermovi scous  dissipation  at  shocks;  these  nonlinear 
phenomena  are  described  using  the  weak-shock  approximation. 
Results  from  the  uncorrected  weak-shock  model,  which  omits  chemical 
relaxation  effects,  are  known  to  deviate  appreciably  from  long- 
range  measurements  [P.  H.  Rogers,  J.  Acoust.  Soc .  Am.  68,  1418-1419 
(  1977)];  the  discrepancy  is  explained  by  showing,  with  the  aid  of 
the  present  model,  that  initial  waveforms  having  different  time 
scales  exhibit  markedly  different  long-range  behaviors.  Specific 
waveform  properties  examined  include  the  peak  pressure,  decay  time 
constant,  positive-phase  impulse,  and  shock  rise  time.  Reasonable 
agreement  is  demonstrated  between  model  predictions  and  long-range 
measurements  in  the  ocean,  as  represented  by  empirical  scaling 
1 aws ." 


★ 

References  cited  in  this  style  indicate  items  in  the  Chronological 
Bibliography,  e.g.,  84-2  means  the  second  entry  in  the  list  for  1984. 


Task  II.  Nonlinear  Prooaqation  in  a  Depth  Dependent  Ocean 


The  effect  of  stratification  of  the  medium  on  finite-amplitude 
distortion  is  considered  in  this  task.  The  approach  is  analytical. 
Nonlinear  distortion  is  basically  determined  by  the  particle  velocity 
amplitude  relative  to  the  local  sound  speed,  the  initial  slope  of  the 
wave,  the  coefficient  of  nonlinearity  8,  and  the  distance  traveled. 
Stratification  greatly  complicates  the  process.  First,  the  amplitude 
depends  on  the  area  of  the  ray  tube  along  which  the  wave  is  propagating, 
and  the  area  in  turn  depends  on  the  acoustical  environment.  The 
changing  environment  also  affects  distortion  directly  through  the  change 
in  sound  speed,  impedance,  and  coefficient  of  nonlinearity.  Despite  the 
complicated  relations  involved,  the  problem  can  be  dealt  with 
analytically.  The  trick  is  to  scale  the  dependent  variable  (pressure) 
and  the  travel  distance  in  such  a  way  that  the  nonlinear  wave  equation 
reduces  to  plane  wave  form.  The  known  analytical  solution  for  the  plane 
wave  problem  may  then  be  used. 

The  primary  assumption  that  makes  the  (relatively  simple) 
analytical  approach  feasible  is  that  the  ray  paths  for  finite-amplitude 
waves  are  the  same  as  for  small-signal  waves.  The  ray  program  used  in 
our  work  is  MEDUSA,  a  very  versatile  routine  developed  by  Foreman.^ 

The  requisite  scalings  are  as  follows  (84-1) 


1/2,  _  ,-1/2 


P  =  Pi 


(reduced  pressure) 


-  f  a/s 

X"  J  ~a\S~)  ds'  (reduced  distance),  (7 

s0  °' 

where  S  is  ray  tube  area,  p  is  density,  s  is  arc  length  along  the  ray 
path,  and  a  is  the  thermodynamic  quantity 

a  =  B(pc5)'I/,?  .  (R 

Subscript  zero  means  a  reference  value,  normally  the  value  at  the  source 
(or  effective  source).  Another  important  quantity  is  the  delay  time 
along  the  path, 


An  important  special  case  is  that  of  a  homogenous  medium,  for 

which,  since  S=4  R2,  Eq.  (7)  reduces  to 

*  =  R0£n(R/R0)  .  (10 

This  is  the  classical  expression  for  the  reduced  distance  for  ordinary 

spherical  waves  of  finite-amplitude  (Task  I).  Because  Morfey  expected 
that  Eq.  (10)  would  give  a  resonably  good  first  estimate  of  Eq.  (7),  he 
introduced  the  quantity  G  by  means  of  the  definition 

x  =  sQ£n(Gs/so)  (11 

(G  may  be  calculated  by  equating  the  right  hand  sides  of  Eqs.  (7)  and 
(11)).  The  departure  of  G  from  unity  provides  a  measure  of  the  extent 
to  which  stratification  affects  distortion. 

Morfey's  technical  report  (84-1)  provides  all  details  of  the 
results  for  Task  II.  Topics  include  conversion  of  Eq.  (7)  to  an 

integral  in  terms  of  ray  coordinates,  exact  analytical  expressions  in 
certain  cases  (for  example,  for  a  linear  sound  speed  profile), 
asymptotic  approximations  for  G,  examples  of  the  calculation  of  G  for 
typical  ocean  profiles,  and  calculation  of  G  when  a  caustic  is 

approached.  Appendices  give  formulas  for  acoustical  properties 
(including  8)  for  the  ocean,  useful  approximate  relations  for  a,  6,  and 
p,  some  asymptotic  ray  relations  for  small  angles,  an  analysis  of 
vertical  propagation,  and  evaluation  of  G  by  means  of  a  ray  path 
integration. 

C.  Merger  of  Tasks  I  and  II.  Applications 


Near  the  end  of  the  contract,  work  was  begun  to  merge  the  results 
of  Tasks  I  and  II,  in  particular,  to  generalize  the  Task  I  computer 
program  to  provide  for  propagation  through  an  inhomogeneous  ocean.  Ray 
paths  involving  reflections  and  caustics  were  avoided.  Although  not 
completed  under  the  present  contract,  the  merger  was  continued  and 
finished  under  Contracts  N00014-75-C-0867  and  N00014-84-K-0574.  The 


computer  program  was  used  to  run  a  number  of  examples.  The  importance 
of  the  various  physical  mechanisms  at  work  --  finite-amplitude 
distortion  by  itself  (both  with  and  without  stratification),  absorption 
and  dispersion  alone,  finite-amplitude  distortion  combined  with 
absorption  and  dispersion,  and  so  on  --  was  investigated.  The  general 
question  that  gave  rise  to  the  research  in  the  first  place,  "Are 
nonlinear  effects  important  in  long  range  propagation?”  was  also 
explored.  A  preliminary,  oral  report  of  this  work  (84-4)  was  followed 
by  a  thesis  and  technical  report  (85-3)  and  subsequently  by  a  paper  at 
an  international  symposium  (86-1).  The  following  excerpt  from  the 
abstract  of  85-3  summarizes  the  results. 

"In  this  report  the  propagation  of  finite  amplitude  acoustic 
signals  through  an  inhomogeneous  ocean  Is  investigated  both 
analytically  and  numerically.  The  effects  of  reflections  and 
focusing  are  not  considered. ..  .In  the  numerical  study  the  ocean  is 
assumed  to  be  stratified.  The  effects  of  inhomogeneity,  ordinary 
attenuation  and  dispersion,  and  nonlinear  propagation  are 
investigated  using  a  numerical  implementation  of  nonlinear 
geometrical  acoustics.  Two  explosion  waveforms  are  considered:  a 
weak  shock  with  an  exponentially  decaying  tail  and  a  more  realistic 
waveform  that  includes  the  first  bubble  pulse.  Numerical 
propagation  of  the  simpler  wave  along  a  58.1  km  path  starting  at  a 
depth  of  300  m  leads  to  the  following  conclusions:  (1)  The  effect 
of  inhomogeneity  on  nonlinear  distortion  is  small.  (2)  Dispersion 
plays  an  important  role  in  determining  the  arrival  time  of  the 
pulse.  (3)  Neither  nonlinearity  nor  ordinary  attenuation  (and 
dispersion)  are  paramount;  both  need  to  be  included.  For  the  more 
realistic  wave  the  propagation  is  along  a  23  km  ray  path  starting 
from  a  depth  of  4300  m.  Two  charge  weights,  0.818  kg  and  22.7  kg 
TNT,  are  assumed.  In  each  case  the  energy  spectrum  of  the  signal 
obtained  by  considering  finite  amplitude  effects  for  the  entire 
23  km  path  is  compared  with  spectra  obtained  by  neglecting  finite 
amplitude  effects  (1)  entirely,  (2)  after  the  first  150  m,  and 
(3)  after  the  first  100m.  Finite  amplitude  effects  are  found  to 
be  of  small  consequence  in  the  case  of  the  0.818  kg  TNT  explosion 


for  frequencies  below  6  kHz  at  distances  beyond  1100  m.  For  the 
22.7  kg  explosion  the  corresponding  quantities  are  4  kHz  and 
1100  m." 


D.  Task  III.  Nonlinear  Propagation  in  a  Caustic  Region 

The  approach  of  a  ray  to  a  caustic  indicates  the  impending 
breakdown  of  the  ray  description  of  the  sound  field.  At  the  caustic 
itself  the  ray  tube  area  vanishes  and  the  pressure  is  predicted  to 
become  infinite.  What  actually  happens  is  that  diffraction,  which  is 
neglected  in  ray  theory,  limits  the  growth  in  amplitude.  Also  not 
included  in  ray  theory  is  the  phase  shift  that  occurs  when  the  wave 
passes  through  the  caustic.  For  small  signals  the  phase  shift  is  tt  for 
a  three-dimensional  focus,  and  ir/2  for  a  two-dimensional  focus  (the 
usual  caustic).  These  results,  if  applied  without  change  to  the 

finite-amplitude  problem,  imply  that  one  need  only  add  a  ff/2  phase  shift 
to  all  frequency  components  of  the  finite-amplitude  wave  when  it  reaches 
the  caustic  (more  accurately,  the  caustic  produces  a  Hilbert  transform 
of  the  incident  wave).  Indeed,  just  such  a  procedure  has  been  used  by 

1C  1C 

Ostrovskii,  Pelinovskii,  and  Fridman,  and  by  Rogers  and  Gardner. 

Buckley  undertook  the  development  of  a  more  accurate  procedure  for 
handling  propagation  through  the  caustic  region.  His  work  is  summarized 
by  the  following  abstract. 

"Nonlinear  acoustic  ray  tracing  in  the  ocean  fails  whenever  a  wave 
front  converges  to  a  caustic  or  focus  and  neighboring  rays 
intersect.  It  is  the  purpose  of  this  paper  to  construct  a 
numerical  algorithm  for  dealing  with  propagation  of  large  amplitude 
pulses  in  the  neighborhood  of  a  caustic.  Introduction  of  a 
time-straining  transformation  of  the  nonlinear  wave  equation 
allows  one  to  attempt  a  perturbation  expansion  of  the  pressure  in 
powers  of  signal  strength.  An  essential  ingredient  of  the  procedure 
is  a  knowledge  of  the  solution  of  the  corresponding  linear  problem. 
An  efficient  procedure  for  calculating  the  linear  solution  for 


incident  N  wave  pulses  having  an  arbitrary  two-dimensional  wave 
front  is  described.  The  corresponding  calculation  of  the  non¬ 
linear  time  straining  is  also  presented.  The  procedure  may  readily 
be  adapted  to  pulses  of  more  general  shape." 

This  abstract  introduced  an  early  version  of  a  manuscript  intended  for 
publication.  Current  plans  call  for  dividing  the  manuscript  in  two. 
The  linear  theory  will  be  the  subject  of  one  manuscript  (86-2),  and 
nonlinear  theory  the  subject  of  the  other  (86-3). 

E.  Nonlinear  Mixing  of  Two  Sound  Beams 

Westervelt  carried  out  an  analysis  of  the  scattering  of  sound  by 
sound.  He  has  long  maintained17  that  interaction  of  two  primary  sound 
beams  produces  no  scattered  difference  frequency  (or  sum  frequency) 
sound  outside  the  interaction  region  unless  the  primary  beams  are 
collinear,  that  is,  unless  the  angle  of  interaction  is  zero.  A  defense 
of  this  position  with  new  calculations  was  presented  in  a  paper  given  at 
an  International  symposium  (84-3). 
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